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(57) -Methode pour construire un module (hydrody- 
namique ou thermodynamique par exemple) desti- 
ne a simuler en temps reel le mode d'ecoulement 
en tout point d'une conduite, d'une veine de fiuide 
polyphasique comprenant au moins une phase li- 
quide et au moins une phase gazeuse, de facon 
qu'il soit le mieux adapte a des conditions operatoi- 
res fixees portant sur un certain nombre de para- 
metres structured et physiques definis relatifs a la 
conduite, et sur un ensemble de grandeurs physi- 
ques definies (grandeurs hydrodynamiques ou 
thermodynamiques par exemple), avec des gam- 
mes de variation fixees pour les parametres et les 
grandeurs physiques. 
- Elle comporte ['utilisation d'un systeme de modeli- 
sation constitue a base de reseaux neuronaux non 
lineaires avec chacun des entrees pour des para- 
metres de structure et des grandeurs physiques, et 
des sorties ou sont disponibles des grandeurs ne- 
cessaires a I'estimation du mode d'ecoulement, et 
au moins une couche interm6diaire. Les reseaux 
neuronaux sont determines iterativement pour 
s'ajuster aux valeurs d'une base d'apprentissage 
avec des tables predefines reliant differentes va- 
leurs obtenues pour les donnees de sortie aux va- 
leurs correspondantes des donnees d'entree. On 
utilise une base d'apprentissage adaptee aux con- 
ditions operatoires imposees et on genere des re- 
seaux neuronaux optimises s'ajustant au mieux aux 
conditions operatoires jmposees. 
- Applications a la moderation d'ecoulements d'hy- 
drocarbures dans des conduites, par exemple. 
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Description 

rOOOIl La presente invention a pour objet une metho- 
de pourformerun module a reseauxneuronaux des m6 

a simuler en temps reel le mode d'ecoulement en tout 
comprenantaumoinsunephaseliquideetaumo,nsune 
chase qazeuse, qui soft le mieux adapte a des cond - 
SoS opStoires et a un ensemble de grandeurs phys, 
ques fixees 
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Etat de la technique 

[0002] L'acheminement des hydrocarbures depuis 
es sftes de production jusqu'aux unites de traftement 
constftue un maillon important de la chain. P^e 
C'est un maillon delicat en raison de la complexity des 
interactions entre les phases constant 
transports Les operateurs ont pour object* prem,e 
d'atteindre une productivity optimaie dans les me. eu 
res conditions de securite. lis dovent done g6rer au 
mieux la Vitesse et la temperature, P ourev «^ d ^; 
tes de charges superflues, des dep&ts mdes.rables et 
des ,rregula 9 rites d'ecoulement. La methode genera.e- 
mentutiHseeconsisteamodeliseraum,euxletransport 

de flux polyphasiques complexes de facon a fournir a 

ferentes parties de la chaTne de product.on, tenant 
compte de la constitution precise de I'effluent, les debts 
etpressionset les modes d'ecoulement. 
[00031 II existe actuellement differents modules log.- 
ciels de simulation du transport de flux polyphas.ques 
complexes, permettant a une stade prteoce de conce- 
voir des equipements de production adaptes 
TO0041 les brevets US 5 550 761 , FR 2.756.044 
US 6 028 992) et FR 2 756 045 (US 5 960 1 87) du de- 
manded, notamment. on connait des modules de mo- 
delisation formant le modele Tacite connu des .gens 
metier, permettant de simuler le transport de flux poly- 
phasiques complexes en regime permanent ou trans,- 
toire et capables de prendre en compte des phenome- 

irrequliere du terrain ou passe la conduite ou de la to- 
po? r aph!edecel.e-ci,que.esspeoi^ 
« terrain slugging » ou « severe slugging » . 
t£Sn La complexite des modules de slmu ation est 
I Primage de celle des phenomenes modelises^ Prec, 
sion e^ performances ne peuvent etre obtenues 
qu'apres un temps de modelisation relativement impor- 
tant qui s'avere difficilement compatible avec une ges- 
tion en temps reel. C'est pourquoi ces modules , de mo- 
delisation ne peuvent en I'etat serv.r a la gestton en 
temps reel de la production. Dans cette opt.que, ,. est 
ainsi apparu necessaire de mettre en oeuvre des me- 
thodes de modelisation permettant un bon comprom.s 
entre rapidite de calcul et precision des resuftats^ 
moos Parlademande debrevetFR 00/09 .889dude- 
mandeur, on connaTt une methode pour eshmer en 



temps reel le mode d'ecoulement en tout point d une 
conduite de structure que I'on peutdefinirpar un certa n 
nombre de parametres de structure, d'une ve.ne de flu - 
de polyphasique definie par plusieurs grandeurs phys.- 
5 ques et ccmprenant des phases liquides et ^gazeuses^ 
Suivant cette methode, on realise une moderation du 
mode d'ecoulement en : 

. en formant un reseau neuronal non lineaire avec 
10 une couche d'entree avec autant d'entrees que de 
parametres de structure et de grandeurs physi- 
ques, une couche de sortie avec autant de sort.es 
que de grandeurs necessaires a I'estimation du mo- 
de d'ecoulement et au moins. une couche 
interm6diaire ; 

en constituant une base d'apprentissage avec des 
tables predefinies reliant differentes valeurs obte- 
nues pour les donnees de sortie aux valeurs cor- 
respondantes des donnees d'entree ; et 

en determinant par iterations des facteurs de pon- 
deration de laf onction d'activation permettant de re- 
lier correctement les valeurs dans les tables des 
donnees d'entr6e et de sortie. 



[0007] De preference, on analyse des donn6es »de 
Lie des neurones de facon a trier, parm. les valeurs 
des donnees de sortie du reseau de neurones, les seu- 
les donnees pertinentes a prendre en compte dans la 
SteZation iterative des coefficients de P onderat.on 
de la fonction d'activation. 

La methode selon I'invention 
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La methode selon I'invention permet de constru.re 
un module (hydrodynamique ou thermodYnamique 
parexemple) destine a simuler en temps reel le mo- 
de d'ecoulement en tout point d'une conduite, d une 
veine de fluide polyphasique ccmprenant au mo.ns 
une phase liquide et au moins une phase gazeuse 
de facon qu'il soit le mieux adapts a des conditions 
operators fixees portent sur un certain n«J*»* 
parametres structures et physiques definis relattfs 
I la conduite, et sur un ensemble de grandeurs phy- 
siques definies (grandeurs hydrodynamiques ou 
thermodynamiques par exemple), avec des gam- 
mes de variation fixees pour les parametres et les 
grandeurs physiques. 

Elle comporte ('utilisation d'un syst&me de modeli- 
sation constitue a base de reseaux neuronaux non 
Iin6aires avec chacun des entrees pour des para- 
metres de structure et des grandeurs phys.ques, et 
des sorties ou sont disponibles des grandeurs ne- 
cessaires a I'estimation du mode d'ecoulement, et 
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au moins une couche intermediaire. Les reseaux 
neuronaux sont determines iterat'rvement pour 
s'ajuster aux valeurs d'une base d'apprentissage 
avec des tables predefinies reliant drfferentes va- 
leurs obtenues pour les donnees de sortie aux va- 
leurs correspondantes des donnees d'entree. 

La methode se caracterise en ce que Ton forme une 
base d'apprentissage adaptee aux conditions ope- 
ratoires impos6es et on g6n6redes reseaux neuro- 
naux optimises s'ajustant au mieux aux conditions 
operatoires imposees. 

[0009] Dans le cas, par exemple, ou le module doit 
etre integre a un module general a la fois hydrodynami- 
que et theimodynamique de simulation d'ecoulements 
polyphasiques, on utilise le modele pour former la base 
d'apprentissage, de manfere a seiectionner i'ensemble 
de grandeurs physiques le mieux adapts au fonctionne- 
ment du module, ainsi que les gammes de variation 
fixees pour les dits parametres et les dites grandeurs 
physiques, et Ton gen&re les reseaux neuronaux opti- 
mises s'ajustant au mieux a la base d'apprentissage for- 
m6e. 

DESCRIPTION DETAILLEE 

Generalites sur la formation du systeme de reseaux 

[0010] On considere une circulation de fluides poly- 
phasiques dans une conduite avec au moins une phase 
liquide et au moins une phase gazeuse et Ton cherche 
a construire un systeme neuronal permettant, a partir 
d'un certain nombre de donnees d'entree g6om6triques 
et physiques relatives a la conduite et de donn6es phy- 
siques sur ies fluides, de donner a chaque instant et 
pour chaque section de la veine fiuide une estimation 
du mode d'ecoulement, ou bien, a partir d'un certain 
nombre de donnees liees aux conditions de pression et 
de temperature sous lesquelles se trouve le fiuide, et a 
la composition du fiuide, de donner le comportement 
thermodynamique du fiuide. L'ensemble de ces don- 
nees d'entree constitue une base de donnees. 

1) entrees et sorties 

[0011] Les donnees d'entree sont par exemple : 

des donnees geometriques et physiques des con- 
duites: diametre, rugosite et angle d'inclinaison du 
conduit, et des donn£es thermodynamiques : mas- 
se voiumique de gaz, masse volumique de liquide, 
viscosity du gaz, viscosity du liquide, etc. ; l'ensem- 
ble de ces donnees constituant un n-uplet ; et 

des donnees hydrodynamiques caracterisant le 
melange: tension superficielle gaz/liquide, fraction 
volumique de gaz, vitesse barycentrique du melan- 



ge, etc. 
ou bien : 

des donnees thermodynamiques relatives a I'envi- 
ronnement du fiuide : pression, temperature, l'en- 
semble de ces donnees constituant un n-uplet, ainsi 
que des donnees caracterisant la composition du 
melange : masse molaire, proprtetes critiques, 
coefficients d'enthalpie, etc. 

[0012] Le module produit par exemple en sortie le 
comportement hydrodynamique des effluents, et, no- 
tamment, le regime d'ecoulement. II evalue et delivre 
sur deux sorties principales des donnees hydrodynami- 
ques dans la partie de conduite dont on souhaite deter- 
miner le type d'ecoulement, la difference dV de Vitesse 
entre gaz et liquide par exemple, ou la fraction p(P £ 
[0 ; 1]) d'ecoulement de type stratifie (type ou la phase 
liquide s'ecoule dans la partie inferieure de la conduite). 
D'autres grandeurs qualifiant le type d'ecoulement peu- 
vent etre calcuiees a partir de ces deux sorties. 
[0013] Dans les applications a la thermodynamique, 
le modele produit en sortie les grandeurs thermodyna- 
miques des effluents, comme, par exemple, le nombre 
et la nature des phases en presence ou les fractions 
molaires des phases. 

2) Structure du systeme de reseaux 

[0014] La fonction que doit modeliser le mod<ele hy- 
drodynamique ou thermodynamique fait appel a diffe- 
rentes sous-fonctions que Ton peut resoudreen utilisant 
si necessaire plusieurs reseaux dont le role au sein du 
module est bien defini. Par exemple, un reseau peut etre 
d6die a la regression de fonction non-lin6aire, alors 
qu'un autre reseau sera dedie a la classification de dif- 
ferents regimes d'ecoulements. Par ailleurs, les liens a 
etablir entre r6seaux ou entre le module et son environ- 
nement exterieur, peuvent necessiter des traitements 
de donnees compiementaires (normalisation des don- 
nees, denormalisation, etc.). Ainsi, on pariera par la sui- 
te de systeme base sur des reseaux de neurones. 
[0015] Pour relier entre elles toutes les donnees hy- 
drodynamiques ou thermodynamiques d'entree et de 
sortie, on construit un systeme de reseaux de neurones, 
ces reseaux etant de preference de type MLP bien con- 
nu des gens de Tart, car il est particuli&rement bien 
adapte a la modeiisation de phenomdnes physiques. Sa 
structure permet de decrire tout a la fois les composan- 
tes dynamiques etstatiques desphenomenes, enfixant 
meme si necessaire certains de ses parametres a une 
valeur reifiee, et done physiquement representative. 
Ainsi, la connaissance d'equations physiques regissant 
les ecoulements permet par exemple d'enrichir le re- 
seau et de I'adapter au mieux aux phenomdnes physi- 
ques qu'il modelise. 

[0016] Le r6seau de neurones comporte (Fig.1) par 
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exemple trois couches : la couche d'entree composee . 
de dix neurones correspondent aux dix donnees (citees 
ci-dessus) du modele physique complet, une couche de 
sortie de deux neurones correspondant aux deux para- 
metres recherches dV et 0, et une couche intermediaire, 
dite couche cachee, dont le nombre de neurones N c est 
optimise. Le reseau est total ement connecte. La non li- 
nearite de ce reseau est obtenue par une fonction d'ac- 
tivation sigmolde regissant le comportement des neu- 
rones de la couche cachee. Les neurones de la couche 
de sortie peuvent etre choisis lineaires ou non-lineaires. 
La fonction d' activation peut etre par exemple la fonction 
identite. 

3) Apprentissage 

[001 7] Les poids de cette structure sont determines a 
Tissue d'une phase d'apprentissage ; au cours de cette 
phase, on nourrit le reseau d'un ensemble de donnees 
constituant la base d'apprentissage du reseau, et on op- 
timise la structure et les poids du reseau en minimisant 
des erreurs constatees pour I'ensemble des echan- 
tillons de la base, entre les donnees de sortie issues du 
calcul du reseau et les donnees attendues a la sortie, 
donnees par la base. Les erreurs peuvent etre les er- 
reurs absolues entre les grandeurs d'entree et de sortie 
ou les erreurs relatives, selon la performance desiree 
pour le reseau. 

[0018] Les facultes de generalisation du reseau sont 
ensuite testees sur sa capacite a bien calculer les deux 
sorties pour des entrees qui lui sont inconnues. 

Particularites de mise en oeuvre 

[0019] Lesysteme base sur des reseaux de neurones 
modelisant le module hydrodynamique ou le module 
thermodynamique est defini par un ensemble de para- 
metres que la phase d'apprentissage permet de deter- 
miner et de fixer. Ces parametres sont adaptables aux 
contraintesfixees par I'utilisation a laquelle est dedi6 ce 
module. 

- Si le module est destine a une utilisation dans un 
cadre general, la base d'apprentissage utilise doit 
etre exhaustive c'est-a-dire Inclure toutes les don- 
nees disponibles concernant les differentes en- 
trees/sorties ; le module est ainsi apte a representer 
I'ensemble des cas que Ton peut rencontrer sur ter- 
rain. 

Si le module est destine a etre utilis6 dans un cadre 
ou un contexte plus specialise, la production d'un 
champ petrolier par exemple, le calage est effectue 
sur une base d'apprentissage plus restreinte con- 
tenant seulement I'ensemble des donnees que Ton 
peut rencontrer sur ce champ uniquement : type ou 
structure des conduites ou pipelines, variations to- 
pographiques, gammes plus restreintes de varia- 



tion de certaines grandeurs hydrodynamiques ou 
thermodynamiques, etc. Les parametres des re- 
seaux neuronaux sont dans un tel cas calcules de 
maniere a obtenir un produit specifique a une con- 
5 figuration donnee, de facon a representer au mieux 
les ecoulements dans des gammes operatoires 
propres a ce champ. 

[0020] Ainsi, ce module est adaptable aux differents 
10 objectifs fixant son utilisation. 

Environnementdu module a reseaux neuronaux 

[0021] Le module developpe peut etre utilise seul, 
is sans autres modules associes. Dans ce cas, il est cale 
sur des bases de donnees constitutes en fonction de 
rutilisation desiree, selon que le cadre est general ou 
specifique. II est independant et permet d'estimer les 
grandeurs hydrodynamiques ou thermodynamiques 
20 susceptibles de permettre le suivi du comportement des 
ecoulements d'effluents dans les pipelines, a partir 
d'une base d'apprentissage regroupant des donnees de 
grandeurs physiques et thermodynamiques fournies 
par ailleurs. 

25 [0022] Le module peut 6ga!ement etre impl6mente 
dans tout modele (ensemble de modules logiciels de si- 
mulation capables de simuler le comportement transi- 
toireet/oustationnairedes ecoulements polyphasiques, 
tel que, par exemple, le modele TACiTE deja cite. Dans 

30 ce tel cas, le module est cale sur des modeles de refe- 
rence representant I'etat de Tart actuel, et sur des bases 
de donnees associees couvrant largement I'ensemble 
du domaine de definition du n-uplet que forment les don- 
nees d'entrees. Son calage tient ainsi compte des spe- 

35 cificites de I'environnement logiciel dans lequel il est 
amene a s'integrer. 

[0023] Pour former la base d'apprentissage la mieux 
adaptee au modele ou le systeme de neurones doit etre 
integre, on utilise avantageusement le modele pour re- 
40 lier entre elles les donnees o°entree et de sortie. 

Exemples d'utilisation 

[0024] Le module a base de reseaux neuronaux con- 
45 vient pourtoute application necessitant une represen- 
tation hydrodynamique d'effluents : realisation decodes 
de simulation predictifs d'ecoulements, de simulation 
temps reel ou de simulations d'entraTn ement, etc. Les 
methodes de controle actif en production peuvent aussi 
so se baser sur un tel module. 

Exemples d'optimisation 

[0025] On considere par exemple le cas ou il s'agit 
55 d'optimiser un reseau donnant la classification du type 
d'ecoulement defluides polyphasiques : ecoulement de 
type stratifie ou de classe C1 (definie plus haut), ecou- 
lement de type intermittent ou de classe C2 (caracterise 
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par une succession de bouchons liquides et gazeux), 
ou bien encore de type disperse ou de classe C3 (le 
liquide etant entraTne sous forme defines gouttelettes), 
en fonction d'entrees comme decrit au paragraphe 
« Generalites sur la formation du reseau ». 
[0026] Dans un premier temps, on forme un reseau 
neuronal a I'aide d'une base generate. La structure op- 
timale resultante que Ton obtient comporte 30 neuro- 
nes, avec une seule couche cachee. La table de confu- 
sion du reseau obtenu (Fig.2) exprime le pourcentage 
de classification obtenu dans les trois classes de sortie 
du reseau pour chaque classe de la base de donnees, 
est mon 

[0027] Dans un second temps, on restreint la base de 
donnees utilisee en apprentissage a quelques configu- 
rations particulieres, qui diminuent les plages de varia- 
tion des differentes entrees du reseau. Le reseau neu- 
ronal optimal est dans ce cas constitue de 28 neurones 
et comporte une seule couche cachee. Le resultat de 
classification est materialise par la table de confusion 
de la Fig.3. 

[0028] On verifie immediatement que les resultats de 
bonne classification, repartis sur la diagonale, sont bien 
ameliores par rapport au cas precedent forme au moyen 
d'une base moins selective et moins adapteo a la con- 
figuration ou le module doit fonctionner. 
[0029] Ainsi, une adaptation de la base de donnees 
aux conditions d'utilisation permet de modifier d'une part 
la structure du reseau, d'autre part d'optimiser les resul- 
tats de generalisation obten us. 
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sees et on genere des reseaux neuronaux optimi- 
ses s'ajustant au mieux aux conditions operatoires 
imposees. 

2. Methode selon la revendication 1 , caracterisee en 
ce que I'ensemble de grandeurs physiques est 
constitue de grandeurs hydrodynamiques. 

3. Methode selon la revendication 1 , caracterisee en 
ce que I'ensemble de grandeurs physiques est 
constitue de grandeurs thermodynamiques. 

4. Methode selon Tune des revendication s preceden- 
tes, caracterisee en ce que, le dit module etant in- 
tegre a un modele general a la fois hydrodynamique 
et thermodynamique .de simulation d'ecoulements 
polyphasiques, on utilise le dit modele pour former 
la base d'apprentissage, de maniere a selectionner 
I'ensemble de grandeurs physiques le mieux adap- 
te au fonctionnement du modele, ainsi que les gam- 
mes de variation fixees pour les dits parametres et 
les dites grandeurs physiques, et Ton genere les re- 
seaux neuronaux optimises s'ajustant au mieux a 
la base d'apprentissage formee. 



Revendications 

1 . Methode pour former un module destine a simuler 35 
en temps reel le mode d'ecoulement en tout point 
d'une conduite, d'une veine de fluide polyphasique 
comprenant au moins une phase liquide et au moins 
une phase gazeuse, de facon qu'il soit le mieux 
adapte a des conditions operatoires fixees portant *o 
sur un certain nombre de parametres structured 
definis relatifs a la conduite, et sur un ensemble de 
grandeurs physiques definies, avec des gammes 
de variation fixees pour les dits parametres et les 
dites grandeurs physiques, au moyen d'un systeme 45 
de moderation constitue a base de reseaux neu- 
ronaux non lineaires avec chacun des entrees pour 
des parametres de structure et des grandeurs phy- 
siques, et des sorties ou sont disponibles des gran- 
deurs necessaires a Pestimation du mode d'ecou- so 
lement, et au moins une couche intermediate, les 
reseaux neuronaux etant determines iterativement 
pour s'ajuster aux valeurs d'une base d'apprentis- 
sage avec des tables predefinies reliant differentes 
valeurs obtenues pour les donnees de sortie aux 55 
valeurs correspondantes des donnees d'entree, ca- 
racterisee en ce que Ton forme une base d'appren- 
tissage adaptee aux conditions operatoires impo- 
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